Describing Sr2Ru04 superconductivity in a generalized Ginzburg— Landau theory 
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We propose a simple explanation of unconventional thermodynamical and magnetic properties ob- 
served for Sr2Ru04. Actually, our two-phase model of superconductivity, based on a straight gen- 
eralization of the Ginzburg-Landau theory, does predict two jumps in the heat capacity as well as 
a double curve for the dependence of the critical temperature on an external magnetic field. Such 
theoretical previsions well agree with the currently available experimental data for Sr2Ru04. 

PACS numbers: 74.20.-z; 74.20.De; 74.70.Pq; 74.81.-g 



In a recent series of papers [l|, 0, H, [3] we have succeeded 
in obtaining a straightforward generalizations of the orig- 
inal Ginzburg-Landau (GL) theory Q in order to de- 
scribe s-wave superconductors endowed with two critical 
temperatures, or even spin-triplet one-phase supercon- 
ductors. The basic idea has been to introduce two differ- 
ent order parameters represented by two charged scalar 
fields (really only one mean field with two distinct gauge 
representations, see below) describing Cooper pairs with 
electrons bound by a weaker or stronger attractive force, 
respectively. The resulting theoretical model is there- 
fore able to describe superconductors with two distinct 
superconducting phases, since the two order parameters 
condensate, in general, at different critical temperatures. 
Peculiar thermal and magnetic properties of these kinds 
of superconductors have been discussed in [3, S] ■ Here we 
only mention that an additional discontinuity in the spe- 
cific heat is predicted, with respect to the conventional 
case, when passing from a superconducting phase to the 
other one. Moreover, at low temperature the London 
penetration length for the superconductors considered is 
strongly reduced, and the Ginzburg-Landau parameter 
K becomes a function of temperature. Instead, in tem- 
perature region between the two phase transition, k is 
constant and the system behaves as a type I or a type 
II superconductors depending on the ratio between the 
two critical temperatures. Such a ratio may be as large 
as v/473 0, [3 (that is, a maximum difference of ^ 15% 
between the two critical temperatures) for very large self- 
interaction of the Cooper pairs with respect to the elec- 
tromagnetic coupling. By allowing a suitable non-linear 



interaction among the two scalar fields, the same theo- 
retical model may as well account for rotational degrees 
of freedom in superconductivity, that is spin-triplet su- 
perconductors with a single phase (the two mutually in- 
teracting order parameters condensate simultaneously at 
a same temperature). In the corresponding model, the 
main thermodynamical and magnetic properties of these 
p-wave superconductors turn out to be essentially the 
same as for the conventional s-wave superconductors. 

All the above-seen properties have prompted us to ex- 
plore the possibility to use the proposed model in or- 
der to understand the intriguing properties exhibited by 
Sr2Ru04 superconductors, which still wait for a com- 
prehensive and solid explanation (see, for example, the 
review in Q). The layered Sr2Ru04 is a superconductor 
with a very low critical temperature and long coherence 
length. It possesses pronounced unconventional features 
such as the invariance of spin susceptibility (across its 
superconducting transition temperature Tc), a strong de- 
pendence of Tc on non-magnetic impurities, evidence for 
two-component order parameter from field distribution 
in the vortex lattice, and so on Q. All these features 
strongly suggest that Sr2Ru04 superconductivity does 
not involve standard s-wave singlet pairing but, rather, 
it could be understood in terms of spin-triplet supercon- 
ductivity. However, the observation of a nodal structure 
of a superconducting gap with a seemingly circular line 
node around a cylindrical Fermi surface, in contrast to 
an expected nodeless p-wave superconductivity in anal- 
ogy to the case of superfluid '^He, together with other 
NMR measurements, leads to a non definitive statement 
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even about the symmetry of the superconducting order 
parameter. In addition to these features, strong (thermo- 
dynamic) evidence for a second superconducting phase 
in Sr2Ru04 exists for high magnetic field appUed along 
a direction H \\ 100. The formation of this second su- 
perconducting state seems closely related to the upper 
critical field limit in a parallel field configuration, as is 
seen in the temperature dependence of Hc2{T) and in 
the scaling breaking of the field dependence of thermal 
conductivity with Hc^ [1, 0] ■ 

Despite the failure in applying standard theoretical 
models to describe such non conventional features of 
Sr2Ru04 superconductivity, the appearance of a second 
discontinuity in the specific heat measurements seems to 
suggest the possible application of our above-mentioned 
GL-like model with two charged scalar fields 0, H, 
Note that, in a Ginzburg-Landau framework, it is quite 
easy to introduce two or more ad-hoc order parameters 
describing different phase transitions, just by considering 
two or more scalar fields with different masses and self- 
interaction coupling constants. In such a way, however, 
the resulting model is not as predictive as the original GL 
theory, since additional unknown physical parameters ap- 
pear. In our present model, instead, the introduction of 
additional degrees of freedom, described by two complex 
scalar fields (/j^, (rather than a single one), does not 
involve new unknown physical constants, so that it is fully 
predictive. As a matter of fact we do require that both 
scalar fields be endowed with equal bare masses m and 
self- interaction coupling constants A. The Lagrangian de- 
scribing the system is, indeed, the following (/i = c = 1) 

C = (i?^(/)w)^^^0w)+m20t0w-^(0Uw)' 

- \f,uF^''- (1) 

where > 0, F^^y = df_,A,j — dyA^ is the electro- 

magnetic field strength, and = 9^-1- 2ieA^ is the 
covariant derivative (2e is the electric charge of a Cooper 
pair). Different critical temperatures Ti and T2 are due 
to different condensations of electrons in Cooper pairs, 
mediated by different effective self-interaction (induced 
by loop renormalization) 0. By starting with the nor- 
mal state system and lowering its temperature, we meet 
a first spontaneous symmetry breaking (SSB) at a criti- 
cal temperature Ti and the medium becomes supercon- 
ducting (phase-I). By further lowering the temperature 
at T = T2, the condensation involving the other order 
parameter is energetically favored and a new phase tran- 
sition starts (phase-II). In a microscopic BCS-like frame- 
work, this corresponds to assume that, as usual, at tem- 
perature Ti electrons at the Fermi surface (with energy 
ep) become bound into Cooper pairs, and the BCS state 
|'0BCs) describing this condensate of electron pairs above 



the filled electron Fermi sea is expressed as follows [8|: 

I^BCS) = W (ww,k Cw -ki + Ww,k ct.kt) |0) > (2) 
k 

where c^i(.|, Cw.-kx are the creation operators of elec- 
tron and hole states around (below) temperature Ti with 
given momentum k, — k (single particle energy Ck) and 
spin up, down, respectively. In the standard BCS ap- 
proximation, the Hamiltonian of such a system is given 
by 

'^ = ^{'^k- ep) ct,kaCw.ka- 
k,CT 

~ Iffeff P ^ Cw,kT^w,-kiCw,-k'iCw,k'T i (3) 
k,k' 

where is the effective clcctron-phonon interaction 
coupling constant at temperature Ti. In the mean field 
approximation, by assuming that each Cooper pair is 
much larger than the typical spacing between particles, 
we can introduce the expectation value 

Aw = -Iffe'ff P $](Cw,-k'iCw,k'T) , (4) 
k' 

and the BCS Hamiltonian becomes approximately 
'^ = ^i^k- ep) ct,kCTCw.k<T + 

k,<T 

+ (cw.kT^l.-kiAw + A^Cw.-kiCw.kt) • (5) 
k 

This Hamiltonian may be easily diagonalized by a Bo- 
goliubov transformation with the introduction of creation 
and annihilation operators b\j, 6w that are quantum su- 
perpositions of electron and hole operators [8| , obtaining 

^ = X] {FvjM ^t.kT^w.kT - -E^w,fc ^'w,-ki^|v,-ki) ' (6) 

k 

the energy eigenvalues being given by k = 
V(efe-ep)2 + |AwP. 

At a lower temperature T2, some of the electrons at 
the Fermi surface already bounded in the formed Cooper 
pairs rearrange themselves to form new Cooper pairs with 
a different effective electron-phonon interaction ij^jj (see 
below), so that below T2 the BCS state is given by 

IV'BCS) = J]^ (uw,kCw,-ki +l'w,kc!|^^kT + 
k 

+Ms,kCs -ki +1's,k4,kT) 1'^) ' 
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with an obvious meaning for the quantities involved. This 
means that an additional gap Ag arises, 



(8) 



and the full Hamiltonian describing the system at tem- 
perature below T2 is then 



w, — kj. 



A peculiar prediction of this model is that two disconti- 
nuities in the specific heat are expected at temperature 
Ti and T2 01, the first one being due to the formation 
of Cooper pairs of the first kind, while the second one 
corresponding to the rearrangement of some electrons in 
Cooper pairs of the second kind. However, in contrast 
with other models with two independent order param- 
eters such critical temperatures may differ at most of 
^ 15% as clear from Eqs. ([T2l) 



,fc^s,-kifol,_kx 



1 < (Ti/Ts)^ < 4/3, 



(13) 



(9) 



where -Bsfe = \/ [ek — ep)^ + |Asp are the energy eigen- 
values for the novel kind of electron-hole states. 

In the most transparent (for our purposes) GL for- 
malism, after the condensations due to the U(l) SSB, 
the mean total free energy results as the sum of con- 
tributions from normal-conducting electrons, and from 
weakly-coupled and strongly-coupled Cooper pairs: 

F = Fn for r>ri, 

F = F^ + a^{T)r,l + ^4 for < T < , 

4 (10) 

F = F^ + a„iTWa + ^Vo + 

+as{T)xo^ + ^xt for T < T2 , 

{riQ is the nonzero expectation value of |(/)wh while xo is 
the nonzero expectation value of Re{(/)s}), with 



aw(T) 



-m 



2 , ^±i£!r2 
16 



as(T) 



-m 



2 , A + 3e2 2 
12 



(11) 

Note that, despite the fact that the bare masses and 
self-interaction coupling constants are the same for both 
scalar fields, the GL effective parameters a^{T), as{T) 
are not: this depending on a different choice for the 
degrees of freedom, described through the scalar fields 
(/>w, 4>8, which undergone condensation P, 0, 0|. In the 
BCS formalism this corresponds to the fact that, at lower 
temperatures, the electron-phonon effective interaction 
may be different than at higher temperatures, due to 
particular constitutive properties of the material consid- 
ered. In fact, even though the bare value of the electron- 
phonon coupling constant 170 does not change with tem- 
perature, the effective coupling g^ff entering the BCS 
Hamiltonian does exhibit such a behavior when intro- 
ducing temperature-dependent radiative corrections, the 
size of the effect depending on the particular system con- 
sidered. The net result is the appearance of two criti- 
cal temperatures, which are defined by the vanishing of 
a^{T) and as{T), that is 



2 ^ 16™^ 
1 A + 4e2' 



r- 



2 _ 12m2 
2 ^ A + 3e2 



(12) 



the value of this ratio depending only on the ratio A/4e2 
of coupling constants of the scalar self-interaction and 
the electromagnetic interaction. Note that, in general, 
the arrangement of the Cooper pairs at temperature T2 
does not automatically guarantee a visible second jump 
in the specific heat exactly at T2, since this process can 
in principle be smoothed over some temperature inter- 
val. It is, however, remarkable that two jumps in the 
heat capacity are just observed for Sr2Ru04 in the pres- 
ence of an external magnetic field. In fact with reference, 
for example, to the measurements reported in 0,0], for 
H > 1.2T a second peak in Ce{T) appears at a second 
critical temperature near the upper limit in (|13p . The 
exact values of these transition temperatures depend on 
the magnetic field applied (they decrease with increasing 
magnetic field) and the effect seems to be suppressed at 
very low temperature, but notably the ratio of the criti- 
cal temperatures appears always close to the upper limit 
7473. 

If the present model will be further confirmed by future 
experiments on other peculiar predictions (see [2, S]), ^ 
possible, general interpretation of Sr2Ru04 superconduc- 
tivity is the following. The striking features of this ma- 
terial strongly depend on the dynamics of two different 
Cooper pairs formed in it, probably due to peculiar ge- 
ometric arrangement of atoms in the compound consid- 
ered . The different effective masses of these pairs (see 
Eqs. (fTTj) ) are responsible of the two peaks in the specific 
heat, corresponding to the condensation of electrons in 
the two different Cooper pairs, which differ one another 
from the relative effective bound of the pairs of electrons. 
The value of the ratio of the two critical temperatures 
close to the upper limit in (|13p . then, points out the 
extremely large value of the scalar self-interaction cou- 
pling A with respect to the electromagnetic coupling (see 
Eqs.lini)). 

In the present model it is assumed that such a scenario 
comes out for electrons grouped in s-wave pairs 0, 0]. 
However, all the known phenomenology which has been 
ascribed to p-wave coupling is seemingly not in con- 
tradiction with the last statement. In fact in Rcf. Q we 
have shown that by switching on a suitable non-linear 
interaction potential among two different Cooper pairs 
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{01, 02 arc the phases of the fields (pi, 4>2) 

m,'A2) = -— In -- = —{9,-62) (14) 

the system behaves as a spin-triplet superconductor, with 
thermodynamical and magnetic properties similar to the 
corresponding ones of the s-wave model (H} . In this way, 
our model is able to reconcile different, unconventional 
features of Sr2Ru04 superconductivity. 

Before concluding, we would address also another in- 
teresting feature reported in the literature, that is the 
"splitting" of the Tc{H) curve [T^]. For sufficiently high 
magnetic fields and low temperatures the first part of 
that curve is ordinarily approximated for small fields to 
a straight line in a GL-like framework 0]. However, as 
results from the previous discussion, for Sr2Ru04, it is 
crucial the condensation of large fields @, 0] , where the 
Tc(H) curve can no longer be approximated to a straight 
line. In such a case, higher eigenvalues (labelled by an 
integer number ti > 0) corresponding to Landau level so- 
lutions of the linearized Ginzburg-Landau equation Q 
come into play in the equation determining the depres- 
sion in transition temperature as a function of the mag- 
netic field. Always referring to s-wave model (1), we 
parameterize this equation as follows (for further details 
see d): 



Tcn{H) ^ Tan " a„(2n + l)fioH 



(15) 



where To„ an are two parameters both depending on 
the elementary flux quantum of vortex lines and on the 
type (weaker or stronger bound electrons) of Cooper pair 
condensed, while the integer n labels the possible en- 
ergy solutions of the linearized GL equation. We have 
analyzed the available data on Tc{H)ioT Sr2Ru04 (for 
definiteness, we took the data from [g, 0]) in the differ- 
ent regimes (different n) pointed out by Eq. ([T5]), and the 
best-fit estimate for the parameters are the following (see 
also Figure 1): 



T^=° = 1.50 - 0.5m 

T:'=^ = 2.6 - 1.52H 

4.03-2.61i? 
4.38 - 2.92i? 



nn=2 



for H < 1.05 

for 1.05 <H< 1.35 

for H > 1.35 
for H > 1.30 



(16) 



(Tc and H are given in K and T, respectively). Notice 
that we can observe a sensible splitting only for n = 2 
eigenstates: in ([TC]) the first equation refers to weakly 
coupled Cooper pairs, the second equation to strongly 
coupled Cooper pairs. 

Note that the ratios of the H coefficients in (|T6| are in 
good agreement with the predictions of Eq. (fT5|) : 2.98 
versus 3, 5.12 and 5.72 versus 5forn = l/ n = 
and n = 2 / n = 0, respectively. Though approximate 
(the approximation of the true curve Tc{H) with a set of 



FIG. 1: Critical temperature versus magnetic field for 
Sr2Ru04 (best fit from the experimental data quoted in 0,01) 



straight lines is, of course, only indicative), our analysis 
allows an immediate recognition of the trigger of the dif- 
ferent regimes ruled by different n and gives some inter- 
esting information on the "splitting" of the Tc{H) curve. 
In fact from (|16|) and, even more, from Fig. 1 it is quite 
evident that such a splitting happens during the transi- 
tion from n = 1 to n = 2, this transition starting at lower 
values of the applied magnetic field for weaker-coupled 
Cooper pairs. The a parameter in Eq. (fTS]) is, indeed, 
inversely proportional to the derivatives, with respect to 
temperature, of coefficients a in PTT]) evaluated at the 
respective critical temperatures [1]: 



(17) 



(here, is the mass of the Cooper pairs), so that, from 
Eqs. pT|) we easily get: 



A 

4 



A 



(18) 



From these relations we deduce that > ag, that is the 
slope of the Tc{H) curve for the weak phase is greater 
that the corresponding one for the strong phase. In par- 
ticular, from Eqs. (fT2|) . the ratio of these two slopes is 
given exactly by the ratio of the two critical tempera- 
tures: 



as 



T2 



(19) 



Thus, the measured slopes for n = 2 give independent 
information on Ti and T2. It is quite remarkable that 
the direct estimate of the ratio T1/T2 from specific heat 
measurements @, 0] is in good agreement with the just 
mentioned ratio of the slopes of Tc{H) for n = 2 in (|16p . 



In conclusion, we have given some evidence that the 
unconventional properties of Sr2Ru04 superconductivity 
may be accounted for by a suitable generalization of the 
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GL theory, as given in Refs. [H, d, H, [l] • In such a model 
the two peaks appearing in the specific heat curve at 
sufficiently high magnetic field applied (and at low tem- 
peratures) correspond to the condensation of two differ- 
ent types of Cooper pairs, where electrons are weaker or 
stronger bound probably due to the peculiar geometry of 
Sr2Ru04 crystals, composed by layers orthogonal to the 
magnetic field. The ratio of the two critical temperatures, 
ruled by the ratio of the scalar self-interaction coupling 
(ruled by A) versus the electromagnetic coupling (ruled 
by e^), is in fact close to the upper limit ■\/4/3 of our 
model, meaning that the self- interaction is in Sr2Ru04 
much larger than the electromagnetic interaction. This 
may be directly inferred from the specific heat measure- 
ments, but it is also in very good agreement with an 
independent estimate from the slopes (for n = 2) of the 
"splitted" Tc{H) curves in the scenario envisaged here. 
The weaker bound Cooper pairs start the transition from 
n=lton = 2ata lower value of the applied magnetic 
field and at a higher temperature, as it is likely expected. 

It is also quite interesting that a trough generaliza- 
tion of the proposed model, accounting for a simple non- 
linear interaction potential among Cooper pairs of differ- 
ent kinds, allows to describe a spin-triplet superconduc- 
tor system, thus recovering even all the unconventional 
features of Sr2Ru04, usually interpreted in terms of p- 
wave pairing. 

Although further experimental and theoretical work is 
desirable for a comprehensive understanding of Sr2Ru04 
superconductivity, the results illustrated here seem 
highly promising to achieve such a goal. 
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Usually, in Sr2Ru04 literature, the splitting is referred 
to the He dependence on the critical temperature(s) Tc. 
However, while for the present purpose the two comple- 
mentary views are identical, we prefer to think at the 
physical situation of the variation of the critical temper- 
ature(s) induced by the applied magnetic field. 
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